Gluon-propagator functional form in the Landau gauge in SU(3) lattice QCD: 
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We study the gluon propagator D^,(x) in the Landau gauge in SU(3) lattice QCD at p = 5.7, 
5.8, and 6.0 at the quenched level. The effective gluon mass is estimated as 400 ~ 600MeV for 
r = (ia) 1 ' 2 = 0.5 ~ 1.0 fm. Through the functional-form analysis of D^(x) obtained in lattice 
QCD, we find that the Landau-gauge gluon propagator D2S(r) is well described by the Yukawa- 
type function e~ mr jr with m ~ 600MeV for r = 0.1 ~ 1.0 fm in the four-dimensional Euclidean 
space-time. In the momentum space, the gluon propagator D^ip 2 ) with (p 2 ) 1 ^ 2 = 0.5 ~ 3 GeV is 
found to be well approximated with a new-type propagator of (p 2 + m 2 )~ 3//2 , which corresponds to 
the four-dimensional Yukawa-type propagator. Associated with the Yukawa-type gluon propagator, 
we derive analytical expressions for the zero-spatial- momentum propagator Do(t), the effective mass 
A/off (i), and the spectral function p(ui) of the gluon field. The mass parameter m turns out to be 
the effective gluon mass in the infrared region of ~ lfm. As a remarkable fact, the obtained gluon 
spectral function p{u>) is almost negative-definite for u) > m, except for a positive 5-functional peak 
at uj — m. 

PACS numbers: 12.38.Aw,12.38.Gc,14.70.Dj 



I. INTRODUCTION 

Quantum chromodynamics (QCD) and the gluon field 
were first proposed by Nambu in 1966 JM just after the 
introduction of color degrees of freedom [2|, and QCD has 
been established as the fundamental gauge theory of the 
strong interaction, through the explanation of asymptotic 
freedom the success of perturbative QCD for high- 
energy phenomena based on theparton model 0, [|[ , and 
lattice QCD calculations d, 0, Hj]. In spite of many suc- 
cesses of QCD, there are still unsolved problems in the 
low-energy region of QCD, owing to its strong-coupling 
nature. Indeed, the mechanism of color confinement is re- 
garded as one of the most difficult important problems in 
theoretical physics 0] , and spontaneous chiral-symmetry 
breaking [lfj is also a difficult issue in describing it quan- 
titatively directly from QCD [111 ]. 

One of the difficulties of QCD lies on the large gap 
between the fundamental fields (quarks and gluons) and 
the observable particles (hadrons). In fact, quarks and 
gluons, which are the building blocks of QCD, have no 
physical asymptotic states, as a result of color confine- 
ment or the gauge invariance. This is in contrast to the 
ordinary perturbation theory, where all the observable 
phenomena can be directly described with the fundamen- 
tal fields appearing in the Lagrangian. Furthermore, the 
field-theoretical description of the confined particles is an 
interesting but unsolved difficult subject. In general, the 
Green function is one of the most basic quantities to de- 
scribe the motions and the interactions of particles [12] , 
but most Green's functions of quarks and gluons are still 



unknown in the nonperturbative description. 

As for the essence of nonperturbative QCD, the central 
issue is gluon dynamics rather than quarks [13j]. In fact, 
the strong gluon interaction makes the QCD vacuum 
highly nontrivial, and color confinement and chiral sym- 
metry breaking are realized even at the quenched level 0] . 
Then, the analysis of gluon properties is the key point to 
clarify the nonperturbative aspects of QCD. In particu- 
lar, the gluon propagator, i.e., the two-point Green func- 
tion is one of the most basic quantities in QCD, and has 
been investigated with much interest in various gauges, 
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such as the Landau gauge 

m m, m m m, p 

the Coulomb gauge [3f 

(MA) gauge [13, [SI] , in the context of various aspects of 
QCD. For example, an infrared-vanishing gluon propaga- 
tor [l5|, [13] is proposed from the mathematical analysis 
of the Gribov horizon 0, [Hj], and the infrared singu- 
larity of the gluon propagator has been investigated in 
terms of color confinement [39j, |4tJ, |4l| , and also from the 
viewpoint of renormalon cancellation p2l |43| . 

Dynamical gluon-mass generation [l3l ] is also an im- 
portant subject related to the infrared gluon propaga- 
tion. In spite of massless perturbative gluons, the gluon 
field is conjectured to acquire a large effective mass as 
the self-energy through the self-interaction of gluons in a 
nonperturbative manner. For, the glueball states, color- 
singlet bound states of gluons, are considered to be fairly 
massive e.g., about 1.5GeV for the lowest ++ and about 
2GeV for the lowest 2 ++ , as indicated in lattice QCD cal- 
culations d,[3- So far, a large effective gluon mass esti- 
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mated about 0.4 ~ 0.8GeV has been proposed in various 
context of QCD physics: analytical framework based on 
QCD [lj|, lattice QCD calculations 0, SH, H2 ~ 
Pomeron physics [461 ] . glueball phenomenology 
and heavy-quark phenomenology [4^ |. Since the color 
SU(iV c ) symmetry is unbroken in QCD, the effective- 
mass generation of gluons is quite different from the Higgs 
mechanism, which is a standard mass generation in quan- 
tum field theories. In the electro-weak unified theory, the 
weak bosons, and Z^, obtain quite a large mass of 
about lOOGeV, as a result of the Higgs mechanism, i.e., 
spontaneous breaking of the SU(2) x U(l) gauge sym- 
metry. However, the gluon mass generation is not a re- 
sult of spontaneous gauge-symmetry breaking, but stems 
from much more complicated nonperturbative dynamics 
of gluons, and generally depends on the gauge choice. 

The nonperturbative effects originate from the strong- 
coupling infrared region of QCD. Actually, a recent lat- 
tice QCD study clarifies that the relevant energy scale 
for confinement is the infrared gluon component below 
1.5GeV [50J: the string tension is almost unchanged even 
by cutting off high-momentum gluons above 1.5GeV. Re- 
flecting the asymptotic freedom or the running coupling 
a s(M 2 )> QCD exhibits various different features accord- 
ing to the energy scale. Here, we roughly classify three 
scale regions of QCD as "ultraviolet (UV)" , "infrared 
(IR) /intermediate (IM)" , and "deep-infrared (Deep-IR)" 
regions, in terms of the length r = (x^x^) 1 ^ 2 . 

• We define the UV region as r < O.lfm, which corre- 
sponds to the high-energy region more than a few 
GeV's. In this region, perturbative QCD is approx- 
imately applicable to the reaction process, and the 
inter-quark potential is almost Coulomb- like [5l| . 
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them and presented high-precision data with the main 
interests in the Deep-IR behavior by using huge-volume 
lattices [32], [3^, [34j . Nevertheless, the functional form of 
the gluon propagator is unclear still now. 

In this paper, we study the functional form of the gluon 
propagator in the Landau gauge in SU(3) lattice QCD 
Monte Carlo calculations, especially for the IR/IM region 
of r = 0.1 ~ l.Ofm, which is considered to be relevant 
for the quark-hadron physics [HI, H3], and also aim to 
describe nonperturbative gluon properties, based on the 
obtained function form of the gluon propagator. 

The organization of this paper is as follows. In Sec. II, 
we give the formalism of the Landau gauge fixing and the 
gluon propagator both in continuum and in lattice QCD. 
In Sec. Ill, we show the lattice QCD result of the gluon 
propagator in the Landau gauge both in the coordinate 
space and in the momentum space. In Sec. IV, we esti- 
mate the effective gluon mass from the gluon propagator 
and the effective-mass plot in lattice QCD. In Sec.V, we 
investigate the functional form of the gluon propagator 
in the Landau gauge, by analyzing the lattice QCD data. 
We show that the Landau-gauge gluon propagator is well 
described with the Yukawa- type function in the IR/IM 
region. In Sec. VI, as the applications of the Yukawa-type 
gluon propagator, we derive analytic expressions for the 
zero-spatial-momentum propagator, the effective mass, 
and the spectral function of the gluon field. Section VII 
is devoted to summary and discussions. 



II. FORMALISM FOR GLUON PROPAGATOR 
IN LANDAU GAUGE 



• We define the IR/IM region as O.lfm ~ r < lfm, 
which ranges from a few hundred MeV to a few GeV 
in energy. In this region, the system is described 
by quark-gluon degrees of freedom in a nonpertur- 
bative way, which is usually substituted with some 
effective models [5j. 

• We define the Deep-IR region as r > lfm, which is 
low-energy below Aqcd ~ 0.2GeV. In this region, 
the perturbative running coupling a s (p 2 ) diverges, 
and the confinement effect is extremely large, so 
that quark-gluon degrees freedom are hidden and 
the system is described by hadrons. 

So far, the gluon propagator has been studied mainly 
in the Landau gauge both in analytic framework [U 
17, " 
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In the UV region, 
the gluon propagator takes a massless perturbative form, 
1/p 2 , apart from the tensor factor. In the IR/IM and 
the Deep-IR regions, the perturbative approach breaks 
down, and we need a nonperturbative approach such as 
lattice QCD calculations Q. Basic important ideas were 
proposed and investigated in the pioneering early-time 



In this section, we review the formalism of the Landau 
gauge fixing and the gluon propagator in the Euclidean 
space-time. 

The Landau gauge is one of the most popular gauges in 
QCD. As a remarkable feature, the Landau gauge keeps 
Lorentz covariance and global SU(A C ) symmetry. Ow- 
ing to these symmetries and the transverse property, the 
color and Lorentz structure of the gluon propagator is 
uniquely determined. 

In the Euclidean space-time formalism such as lattice 
QCD, the Landau gauge is usually defined so as to min- 
imize the gauge-field fluctuation. We then expect that 
only the minimal fluctuation of the gluon field survives 
in the Landau gauge, and the physical essence of gluon 
properties can be investigated without suffering from 
large stochastic fluctuations of gauge degrees of freedom. 



A. Landau gauge fixing 

To begin with, let us consider the Landau gauge fix- 
ing and its physical meaning in Euclidean QCD. In the 
SU(A C ) continuum QCD, the gluon field is expressed as 
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A^(x) = A^(x)T a £ su{N c ) with the generator T a (a = 
1, 2, . . . , iV c 2 - 1) and A*(x) G K. The Landau gauge is 
usually defined by the local condition on the gauge field 
as 

dpApix) = 0. (1) 

In Euclidean QCD, the Landau gauge has a global defi- 
nition to minimize the global quantity, 

Ree J d^x Tr{A,(x)A^x)} = \J d 4 xA«(x)A«(x), 

by the gauge transformation. This global definition is 
more strict, and the local condition (Q} is derived from 
the minimization of R. The global quantity R can be re- 
garded as the total amount of the gauge-field fluctuation 
in the Euclidean space-time. In the global definition, the 
Landau gauge has a clear physical interpretation that it 
maximally suppresses the artificial gauge- field fluctuation 
relating to the gauge degrees of freedom. 

In lattice QCD, the theory is formulated on the dis- 
cretized space-time 0. The QCD action is constructed 
from the link-variable U^x) £ SU(iV c ), instead of the 
gauge field A fl (x) £ su(N c ). The link- variable is defined 
as Uu{x) = e lagA ' i ^ x \ with the lattice spacing a and the 
gauge coupling constant g. The gauge transformation of 
the link- variable is given by 

t^(as)-»n(a;)l7 M (a;)nt(a; + A), (3) 

with the gauge function Q(x) £ SU(N C ). 

In lattice QCD, the Landau gauge fixing is also ex- 
pressed in terms of the link- variable: the Landau gauge 
is defined by the maximization of 

i?iatt = ^^RcTr[W, ( 4 ) 

x fl 

by the gauge transformation of the link- variable. 
For small a, using the expansion 

U M (x) = 1 + tagA^x) - l -a 2 g 2 Al(x) + 0(a 3 ) (5) 
in terms of lattice spacing a, i?i a tt is expressed as 

Aatt = ~^fj2 A » A » + °(° 4 ) + COnSt - ( 6 ) 

x 

Therefore, the maximization of i?i a tt corresponds to the 
minimization of the gauge-field fluctuation as well as the 
continuum theory. This minimization of the gluon-field 
fluctuation in turn justifies the expansion in Eq.([5]). 

B. Gluon propagator in lattice QCD 

In this subsection, we formulate the gluon propagator 
in the Landau gauge. To begin with, we extract the 



gluon field from the gauge-fixed link- variable, which is 
obtained by the gauge transformation of the link- variable 
to maximize i?i a tt ■ Using the expansion ([5]) , we define the 
bare gluon field A^ are (x) as 

(7) 

where the second term is added to make A^ aro traceless. 
The renormalized gluon field A r ° n (x) is obtained by mul- 
tiplying a real renormalization factor Z 3 x ^ 2 as 

A™(x) ee Z- 1/2 A h ^(x). (8) 

We abbreviate A r ° n (x) as A^{x) hereafter. The gluon 
field Afj,(x) defined above is traceless and hermite, and is 
expressed as A^(x) = A a ^(x)T a £ su(N c ) with A°(x) £ 
R. In the Landau gauge which maximizes -Riattj the gluon 
field A^(x) satisfies the local condition, 

dMx) = o, (9) 

with the forward or backward derivative <9 M on the lattice. 

The gluon propagator Dfy(x) is defined by the two- 
point function as 

Df v {x,y) = (A°(x)AUy)) = Df v {x-y). (10) 

Note that the time-ordered product is unnecessary in the 
Euclidean metric, and the translational invariance of the 
vacuum leads to the (x — y)-dependence. From Eq.©, 
(x — y) satisfies the transverse property, 

dlDf v {x-y) = dy,Df v {x-y)=Q. (11) 

Next, we consider the gluon propagator D°^ v {jp) in the 
momentum space, which is defined by the Fourier trans- 
formation of the coordinate-space propagator as 

D$(p) ee J d±x e**D$(x). (12) 

On the Li x L 2 x L3 x L4 lattice, this Fourier trans- 
formation is discretized, and the momentum-space gluon 
propagator is expressed as 

= ( 13 ) 

X 

Here, the discretized momentum p., and the continuum 
momentum p^ are defined as [1, HE H3] 

27m M 2 . fpu,a\ 2 . / 7m„\ 

P ^a Sin [-T) = ~a Sm {-L;)> 

(14) 

with n M = 0, 1, 2, - 1. 

In the Landau gauge, the color and tensor structure of 
the propagator is uniquely determined as 

£%<p)=Dtf)S*(6r-2j&y (15) 
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from the SU(iV c ) global symmetry, the Lorentz symme- 
try, and the transverse property, 



(16) 



Therefore, we only have to consider the scalar factor of 
the gluon propagator, D{p 2 ), which is a function of the 
continuum- momentum squared p 2 — p a p a - 

In the coordinate space, we investigate the scalar com- 
bination of the gluon propagator 



D(r) = 



1 



-DH{x) = 



1 



-K(z)^(0)>, 
(17) 



3(N 2 - 1) wv ; 3(N 2 - I.) 

as a function of the four-dimensional Euclidean distance, 

r = \x\ = (x^x^) 1/2 . (18) 

Here, the denominator factor 3 in Eq. (|17p has been intro- 
duced considering the tensor factor S^p — p^p^/p 2 = 3. 

The scalar factor D(p 2 ) in Eq. ([T5|) is expressed by the 
Fourier transformation of D(r) as 



D(p 2 ) = 



3(iV c 2 -l)- 



from which one can prove that D (r) depends only on r 
near the continuum limit. In this paper, we call D(r) 
and D(p 2 ) "scalar- type propagator". 



III. LATTICE QCD RESULT FOR GLUON 
PROPAGATOR 

We perform SU(3) lattice QCD Monte Carlo calcula- 
tions at the quenched level using the standard plaquette 
action. Here, we adopt three different lattices with the 
lattice parameter j3 = 2N c /g 2 =5.7, 5.8, and 6.0. The 
used lattice size is 16 3 x 32, 20 3 x 32, and 32 4 at /3 = 5.7, 
5.8, and 6.0, respectively. In this lattice calculation, we 
mainly use (3 = 6.0. 

The lattice spacing a is found to be a = 0.186,0.152, 
and 0.104fm, at (3 = 5.7, 5.8, and 6.0, respectively, when 
the scale is determined so as to reproduce the string ten- 
sion as \fa = 427MeV from the static QQ potential [5~j| . 

The number of used gauge configurations is 50, 40 and 
30 for (3 =5.7, 5.8, and 6.0, respectively. The gauge con- 
figurations are picked up every 1,000 sweeps after a ther- 
malization of 20,000 sweeps. We summarize the param- 
eter and calculation conditions in Table |TJ 

Here, we briefly explain the actual procedure of the 
gluon-propagator calculation. For each gauge configura- 
tion, we perform the Landau gauge fixing by the gauge 
transformation to maximize i?i a t defined in Eq.(j4|), and 
obtain the gluon field A^x) defined in Eq.© from the 
gauge-fixed link- variable. Then, we construct the scalar- 
type gluon propagator D(r) from the two-point function 
of the gluon field, as shown in Eq. (fTT)) . Finally, we calcu- 
late the momentum-space gluon propagator D(p 2 ) using 
the discrete Fourier transformation, as shown in Eq. (|19[) . 



TABLE I: The lattice parameter /3, lattice size, and the gauge- 
configuration number 7V con f. The corresponding lattice spac- 
ing a and the lattice volume in the physical unit are added. 
The lattice spacing a is determined so as to reproduce the 
string tension yfo — 427MeV. 



(3 


Lattice size 


a [fm] 


Volume [fm 4 ] 


N 

1 ' cont 


5.7 


16 3 x 32 


0.186 


2.976 3 x 5.952 


50 


5.8 


20 3 x 32 


0.152 


3.040 3 x 4.864 


40 


6.0 


32 3 x 32 


0.104 


3.328 3 x 3.328 


30 



As for the overall renormalization factor Z3 , the gluon 
field {x) is renormalized so as to make the gluon prop- 
agator D{p 2 ) coincide with the tree- level propagator l/p 2 
at /I [27|, «-e., 



D(p 2 ) 



1 

7- 



(20) 



In the calculation of the gluon propagator, we take 
the advantage of the translational symmetry to improve 
statistics, and we adopt the jackknife method to estimate 
the statistical error. 

Figure Q] and [2] show the lattice QCD results at /3=5.7, 
5.8, and 6.0 for the scalar-type gluon propagator D{r) = 
D™(x)/24: and D{p 2 ) = r)^(p)/24, respectively. These 
figures include not only on-axis data but also off-axis 
data. In these figures, the statistical error is rather small 
in the depicted region, and the statistical error bars are 
hidden in the symbols. Here, we choose the renormaliza- 
tion scale at p, = 4GeV for (3 — 6.0 [271 ]. and make cor- 
responding rescaling for /3=5.7 and 5.8. We thus obtain 
D(r) as a single- valued function of the four-dimensional 
Euclidean distance r = (xaXa) 1 ^ 2 . 




0.4 0.6 
r [fm] 



FIG. 1: Lattice QCD results of the scalar-type gluon prop- 
agator D(r) = ELi El=i <A£(a:)A£(0)>/24 as the function 
of the four- dimensional Euclidean distance r = (xaXa) 1 ^ 2 in 
the Landau gauge at /3 = 5.7, 5.8, and 6.0. 
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FIG. 2: Lattice QCD results of the scalar-type gluon prop- 
agator D(p 2 ) = Y, x e i§ ' x D{r) plotted against p = (p M p M ) 1/2 



- sin - 

a V 



in the Landau gauge 



with the momentum 
at /3 = 5.7, 5.8, and 6.0. We renormalize the propagator to 
satisfy the renormalize condition D(p 2 )\ p 2 =)1 2 = l/p 2 at the 
scale n = 4GeV. The dash-dotted line denotes the tree-level 
massless propagator, 1/p 2 . 



The lattice result of D(p 2 ) is almost a single- valued 
function of the magnitude of the continuum momentum, 
p = (paPa) 1 ^ 2 - We confirm that our lattice QCD result 
of D(p 2 ) is consistent with that obtained in the previous 
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12911 . although recent huge- volume 
[3J] indicate a suppression of the 
gluon propagator in the Deep-IR region (p < 0.5GeV), 
compared with the smaller lattice result. The lattice data 
of D(p 2 ) seem to be consistent with the tree-level mass- 
less propagator D tlee (p 2 ) = 1/p 2 for large p 2 , but their 
behaviors are largely different in the low-energy region 
below a few GeVs. 



form Lagrangian [53[ in the Euclidean metric, 

C = \ faAi - d v A^f + \rr?AlAl - ^ (d^f , 

(21) 

where the parameter a = corresponds to the Landau 
gauge. The propagator of the massive vector field A a is 
derived from the Lagrangian as 



i ^( L Ji + ")w^ (22) 



D ab (v) = 



p* — am* 



Taking a = 0, the massive- vector propagator in the Lan- 
dau gauge is obtained as 



1 



b jrab I x 

2_i™2° „2 



p z + mr 



PlJ.Pu 



(23) 



This propagator satisfies the transverse property of the 
Landau gauge, p^D^ u {p) = p v D ab {p) = 0, which corre- 
sponds to the condition, d^A"(x) = 0. 

In this case, the scalar-type propagator D(p 2 ) reads 

1 



d{p 2 ) = ^n™(p) 



p 2 + m 2 



(24) 



and its Fourier transformation gives the scalar-type prop- 
agator D(r) in the coordinate space as [37j 



D(r) = 



" „— was n 



1 m 



(2tt)4 



47r z r 



where Ki{mr) is the modified Bessel function. The 
derivation of this formula is shown in Appendix A. For 
large r, K\{mr) behaves asymptotically as 



K\ (mr) 



2mr 



(26) 



IV. EFFECTIVE GLUON MASS 

In this section, we investigate the effective gluon mass 
in the Landau gauge using the gluonic correlation ob- 
tained in lattice QCD. We derive the massive- vector 
propagator in coordinate space, and estimate the effec- 
tive gluon mass, by comparing the lattice gluon propa- 
gator with the massive propagator. Also, we investigate 
the effective-mass plot of gluons obtained from the zero- 
spatial- momentum propagator in lattice QCD. 



A. Comparison with massive propagator 

We first derive the free massive- vector propagator form 
in coordinate space as a useful guide to analyze the gluon 
propagator in lattice QCD. Here, we use the Stueckerberg 



and therefore the massive propagator behaves as D(r) ~ 

r -3/2 g -mr 

Now, we compare the scalar- type gluon propagator 
D(r) obtained in lattice QCD with the massive propa- 
gator, and estimate the effective gluon mass through the 
fit-analysis. Considering the functional form of the mas- 
sive propagator as shown in Eq. (|2"5|) , we here adopt the 
fit-function defined by 



D maSB (r) = A—Ki(mr), 
r 



(27) 



with the mass parameter m and a dimensionless param- 
eter A. 

For example, we show in Fig|3] the fit result of the 
lattice data D(r) at f3 = 6.0 with D mass (r) in the fit- 
range of r = 0.6 ~ l.Ofm. In this range, this fit seems 
well, and the effective mass m is estimated to be about 
500MeV from this fit. However, FigJ3] shows that the 
lattice gluon propagator D(r) cannot be described with 
-DmassM in the whole region of r = 0.1 ~ l.Ofm. 
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0.2 0.4 0.6 

r [fm] 



FIG. 3: A typical example of the fit-analysis of the lattice 
gluon propagator D(r) with the fit-function -D maS s(r) of the 
massive-vector propagator denoted by the dashed line, The 
fit is done for the lattice data at p = 6.0 in the fit-range of 
r = 0.6 ~ l.Ofm. 



We investigate the fit-analysis for the lattice gluon 
propagator D(r) with D mliss (r) for several fit-ranges, and 
the results are summarized in Table |TTJ From this fit- 
analysis in the Landau gauge, the effective gluon mass 
is estimated as m = 400 ~ 600MeV in the infrared re- 
gion of r — 0.5 ~ l.Ofm, although there is a significant 
r-dependence of m, i.e., m is small at short distances. 



TABLE II: The fit result for the scalar-type gluon propaga- 
tor D(r) obtained in lattice QCD at /3 = 5.7, 5.8, and 6.0. 
The fit function is the massive- vector propagator D mass (r) = 
Amr -1 Ki(mr) . The best fit parameters (m,A) are listed to- 
gether with the fit range and \ 2 /N&t. The fit data with large 
X 2 /Adf are omitted. 





fit ran{ 


;e [fm] 


m [GeV] 


A 


X 2 /N dC 


6.0 


0.5 - 


' 0.7 


0.405(7) 


0.094(2) 


2.07847 




0.6 - 


' 0.8 


0.475(6) 


0.112(2) 


0.588604 




0.7- 


' 0.9 


0.551(8) 


0.140(4) 


0.454586 




0.8 - 


' 1.0 


0.582(10) 


0.158(6) 


0.340416 




0.6 - 


' 1.0 


0.517(5) 


0.125(2) 


1.67541 


5.8 


0.6 - 


' 0.8 


0.502(10) 


0.118(4) 


1.68858 




0.7- 


- 0.9 


0.549(15) 


0.138(7) 


0.773316 




0.8 - 


- 1.0 


0.576(12) 


0.151(7) 


0.414934 




0.6 - 


' 1.0 


0.526(6) 


0.126(2) 


0.576734 


5.7 


0.7- 


' 0.9 


0.626(48) 


0.162(26) 


3.76344 




0.8 - 


> 1.0 


0.618(32) 


0.159(19) 


2.26721 



B. Effective- mass plot of gluons 

In the previous subsection, we estimate the effec- 
tive gluon mass from the fit-analysis for the lattice 



gluon propagator. Now, we estimate the effective gluon 
mass from the effective-mass analysis with zero-spatial- 
momentum propagator Do(t) in the Landau gauge 
I20I l2ll . I22I I23I ]. This method is often used for hadrons 
as a standard mass measurement in lattice QCD [8|. For 
the simple notation, we use the lattice unit of a = 1 in 
this subsection. 

We define the zero-spatial-momentum propagator 
Do(t) of gluons as 



D (t) 



1 

24 



^2{Al(x,t)A^(6,0)) =J2D(r), (28) 



where the total spatial momentum is projected to be zero. 
Using the translational invariance, Do(t) can be rewritten 
as the wall-to-wall correlator, 



D (t) 



1 



;{5>2(*,t)}{£ A 2Gfc°)}> 



(29) 

In the actual lattice QCD calculation, we adopt the wall- 
to-wall correlator to improve statistics with an easy task. 
The effective mass of gluons is defined by 



M eS (t) = ]n{D (t)/D (t + l)}, 



(30) 



in the case of large temporal lattice size. In the numerical 
analysis, we take account of the temporal periodicity used 
in lattice calculations. On the lattice with the temporal 
size N t (— Li), Do(t) is expected to behave as 

D (t) oc e" mt + e-™^*-*) oc cosh - t)}, (31) 

and therefore we define the effective mass M e g (t) by 
D {t + 1) cosh [M rf (t) (N t /2 - (t + 1))] 



D (t) 



cosh [M eS (t) (JVt/2-t)] 



(32) 



which is reduced to Eq. ({3"0"l) in the large N t limit. 

Figure [4] shows the plot of the effective mass M c g(t) 
of gluons in lattice QCD at j3 = 6.0. The effective mass 
M Q g(t) of gluons is an increasing function for small t, and 
approximately constant for t = 0.4 ~ l.Ofm, and the ef- 
fective gluon mass is estimated to be about 500MeV from 
the value of M e g(t) in the range of t = 0.4 ~ l.Ofm. This 
tendency and the estimated value are consistent with the 
previous results [H, [2(| HH, HI] and those obtained 
from the analysis of the gluon propagator in the previous 
subsection. 

Note here that the effective mass M e ff(t) of glu- 
ons exhibits an increasing behavior, which is unusual 
[H, HI EH, HIH. In the usual hadron-mass calcula- 
tion, the effective mass is always a decreasing function of 
t, due to the positive contribution from physical excited 
states at the short distance. The mathematical proof of 
this statement is as follows: The zero-spatial-momentum 
propagator G(t) is expressed in the Lehmann-Kallen rep- 
resentation as 



G(t) = E« 



(33) 
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FIG. 4: The effective mass M e g(t) of gluons in the Landau 
gauge in lattice QCD at /3 = 6.0, i.e., a = 0.104fm. 



with the spectral weight Cj. In the continuum formalism, 
the effective mass is expressed as M(t) = — ^ lnG(i) and 
satisfies 



eft 2 

')(E< 



(34) 



Cim t e 



(E< 



If all the spectral weights are non-negative as Cj > 0, 
this is always non-positive, i.e., 4rM(t) < 0, due to 
the Cauchy-Schwartz inequality [12], so that the effec- 
tive mass M(t) is generally a decreasing function in all 
the region of t. 

However, the effective mass M e s(t) of the gluon ex- 
hibits an anomalous increasing behavior, and this fact 
in turn indicates that t he g luon s p ectral function is not 
positive-definite [HI, [2(J l2ll. l22l [23j due to its unphysical 
nature. Note also that, from Eq. (j3"4]) . we can formally 



obtain 4iM(t) < 0, if all the spectral weights are non- 
positive as Ci < 0. In fact, the increasing property of 
M(t), i.e., 4iM{t) > 0, can be realized, only when there 
is some suitable coexistence of positive and negative val- 
ues in the spectral weights Cj. 

Here, we summarize the lattice QCD result of the effec- 
tive gluon mass in the Landau gauge. The effective gluon 
mass exhibits a significant scale-dependence, and it takes 
a small value at short distances. Quantitatively, the effec- 
tive gluon mass is estimated to be about 400 ~ 600MeV 
in the infrared region ~ l.Ofm. This value seems consis- 
tent with the gluon mass suggested by Cornwall 



V. FUNCTIONAL FORM OF THE GLUON 
PROPAGATOR IN THE LANDAU GAUGE 

In this section, we study the functional form of the 
duon propagator in the Landau gauge in SU(3) lattice 



QCD. In the high-energy region, we already know the ap- 
plicability of perturbative QCD, where gluons are mass- 
less. The perturbative gluon propagator is simply de- 
scribed with 1/p 2 in the covariant gauge, similar to the 
photon propagator in QED. In the infrared region, how- 
ever, the gluon is expected to acquire a large effective 
mass due to nonperturbative QCD effects, as was also 
indicated in the previous section. In fact, the functional 
form of the gluon propagator can be changed according 
to the scale. 

So far, the functional form of the gluon propagator has 
been studied both in analytical framework [lil, fla. fl7t l39j 
and in lattice QCD [H, [H EE HI . In the UV region, 
the lattice QCD studies have shown the perturbative be- 
havior of the gluon propagator. In the Deep-IR region, 
the gluon propagator and its behavior have been investi- 
gated with the theoretical interest of "infrared vanishing" 
behavior in the context of the Gribov horizon, "dipole 
(1/p )-singularity" in the relation to color confinement, 
and so on. Here, we aim to determine the functional form 
of the gluon propagator in the infrared/intermediate re- 
gion of r = (x a x a y/ 2 = 0.1 ~ l.Ofm, which is the rele- 
vant scale of quark-hadron physics. 



Functional form candidates 



In the previous section, we compare the gluon prop- 
agator with the massive propagator D mass (r), which 
corresponds to (p 2 + m 2 ) -1 in the momentum space. 
However, as shown in Fig[31 it is difficult to reproduce 
with D m!iSS (r) the lattice result in the whole region of 
r = 0.1 ~ l.Ofm. In FigEl we choose D mass (r) appro- 
priate for a large-r region, but it gives too large value 
than the lattice result in the small r region. In fact, at 
the short distance, the gluon propagator shows a larger 
reduction than the simple massive propagator. Then, 
one may think of a larger mass at shorter distance in 
the form of (p 2 + to 2 ) -1 . But, a smaller effective mass 
is obtained at the short distance, as was shown in the 
previous section. This indicates that the functional form 
of the gluon propagator itself is largely changed from the 
simple massive propagator D mass {r). 

As the candidate form of the gluon propagator in the 
momentum space, we here consider (p 2 + to 2 ) -3 / 2 and 



as well as the massive propagator, (p 



(p 2 

m 2 ) 1 . Both of two candidates qualitatively satisfy the 
above-mentioned behavior, i.e., larger reduction at the 
short distance. 



For the actual analysis of the lattice QCD result, 
we mainly consider the gluon propagator in the coor- 
dinate space instead of the momentum space, since the 
coordinate-space variable is more directly obtained in lat- 
tice QCD. 
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1. Yukawa-type propagator 

First, we consider (p 2 + m 2 )~ 3 / 2 type propagator. In 
the coordinate space, this corresponds to the Yukawa- 
type function, since the Fourier transformation in four- 
dimension Euclidean space is given by 

f d 4 v 1 11 

/ _^_J_ ipx * _ - 1 - -ror for\ 

J (2tt) 4 ' (p 2 + to 2 ) 3 / 2 4tt 2 r ' V ' 

We show the derivation of this formula in Appendix A. 

Then, we call this form "Yukawa- type function" or 
"Yukawa-type propagator". Usually, the Yukawa-type 
function is obtained by the three-dimensional Fourier 
transformation of (p 2 + to 2 ) -1 . It is notable that this 
Fourier transformation (I35[) is calculated in the four- 
dimensional Euclidean space-time, and therefore the 
momentum-space function takes an unfamiliar form as 
(p 2 + m 2 )~ 3 / 2 . 

For the analysis of the gluon propagator, we introduce 
the definite form of the Yukawa-type propagator as 

777 

^Yukawa(r) = A-e~ mr , (36) 
r 

with a "mass" parameter m and a dimensionless parame- 
ter A. With this form, we analyze the lattice QCD result 
of the gluon propagator in the coordinate space. 

The Fourier transformation of -DyukawaM is given by 

^Yukawa (p 2 ) = J d 4 xe lpx D Yukawa (r) = 4 ^ ^3/2 ■ 

This functional form will be used for the analysis of the 
gluon propagator in the momentum-space in Sec.V-C. 

2. Dipole-type propagator 

First, we consider (p 2 + m 2 )~ 2 type function, which 
we call "dipole-type" . In the coordinate space, this func- 
tion corresponds to the modified Bessel function K (mr) , 
because of the four-dimensional Fourier transformation, 

/(l^pip^^H- (38) 

The derivation of this formula is shown in Appendix A. 
We define the dipole-type propagator Ddi po ie{' r ) as 

AiipoicM = Am 2 K (mr), (39) 

with a "mass" parameter m and a dimensionless param- 
eter A. From the asymptotic form of the modified Bessel 
function, -DdipoleO") behaves as 

^W^m^le-, (40) 
asymptotically for large r. 



3. Summary of three fit-functions 

Here, we summarize the three fit-functions in Table HTT1 
In the coordinate space, the functional forms look rather 
different, i.e., K\(mr)/r, e~ mr /r, and Ko(mr). How- 
ever, there are systematic relations in their asymptotic 
form and their momentum representation, as shown in 
Table IIII1 In the coordinate space, the difference of the 
asymptotic form is just the power of the prefactor as 
r -n/2 ( n = i ( 2,3). 

TABLE III: Summary of the functional form candidates, 
Amuse (r), -DyukawaM, and Ajipoic (r) , together with their 
asymptotic form and their momentum representation. 



Functional form Asymptotic form Momentum space 
Dmass r-'ftH r- 3/2 e- mr (p 2 + m 2 )' 1 

Oyu^ r^e" mr r~ x e- mr (p 2 + m 2 )" 3 ' 2 

Odipoic K {mr) T -^e- mr (p 2 + m 2 )~ 2 



From the aspect of the space-time dimension, the 
Yukawa-type and dipole-type functions may be regarded 
as "low-dimensional" propagator forms. As was al- 
ready mentioned, the Yukawa-type function has a three- 
dimensional character, since the Yukawa function is ob- 
tained by the Fourier transformation of the massive prop- 
agator (p 2 + m 2 ) -1 in the three-dimensional space-time. 
From this viewpoint, the dipole-type function as Ko(mr) 
has a two-dimensional character, since this function is 
obtained by the Fourier transformation of (p 2 + m 2 )" 1 in 
the two-dimensional space-time as 

f d 2 v 1 1 

/ TTvi^ X ^r—-2 = 7rK,{mr). (41) 
J {2ny p z + vn 2n 

The derivation of this formula is shown in Appendix A. 



B. Comparison of lattice QCD results with 
fit-functions 

In this subsection, we compare the lattice gluon propa- 
gator with the three fit-functions, D masa (r), -Dyukawa (*■■), 
and -DdipoicM- In the analysis, we consider the scalar- 
type gluon propagator D(r) in the Landau gauge in the 
range of r = 0.1 ~ l.Ofm obtained in SU(3) lattice QCD, 
and try to reproduce the lattice data through the fit- 
analysis with various range of r for each fit-function. 

Figure [5] shows the typical example of the fit result of 
the lattice gluon propagator D(r) at (3 = 6.0. For each 
fit-function, the best-fit parameters (m, A) and the fit 
range are listed in Table ITVl 

As a remarkable fact, the Yukawa- type function 
^YukawaM well reproduces the lattice QCD data in the 
whole region of r = 0.1 ~ l.Ofm. On the other hand, 
the dipole-type function Ddipoie(^) fads to reproduce the 
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0.2 0.4 0.6 0.8 1 

r [fm] 

FIG. 5: The fit result of the lattice gluon propagator D(r) 
with the three functional forms, -DYukawa(Y) (solid line), 
AiipoicM (dotted line), and -Dmass(r) (dashed line). The 
Yukawa- type function £>Yukawa(r) well reproduces the lattice 
result at (3 — 6.0 in the whole region of r = 0.1 ~ l.Ofm. 



TABLE IV: The best-fit parameters (m, A) and the fit range 
in the fit analysis of the lattice gluon propagator D(r) at 
(3 = 6.0 with the three functions, D mass (r) = AmK\{mr) /r, 
-DYukawa(r) = Ame~ mr /r, and Aiipoie(r) = Am 2 K (mr). 



Functional form 


Fit 


range [fm] 


m [GeV] 


A 




0.6 > 


- 1.0 


0.517(5) 


0.125(2) 


■^Yukawa 


0.1 r 


- 1.0 


0.624(8) 


0.162(2) 


-^dipole 


0.4' 


- 1.0 


0.817(1) 


0.123(1) 



whole region of the lattice data, since it gives too strong 
reduction at the short distance. The massive propagator 
-Dmass(^) also fails to reproduce the whole region of the 
lattice data, as was already shown. 

To see the difference of the three fit results clearer, we 
show in FigJS]the ratio of the lattice QCD data £>iatt (r) to 
the three fit-functions on the scalar-type gluon propaga- 
tor, i.e., Aatt/Anass, Aatt / ^Yukawa , and -Dlatt/Aiipole- 

One finds Matt /^Yukawa — 1 in the whole region of 
r = 0.1 ~ lfm, while Aatt/Anass and Aatt/Alipole differ 
from the unity for small r. 



C. Yukawa- type gluon propagator in Landau gauge 

We thus find an appropriate functional form of the 
gluon propagator in the Landau gauge. In fact, the 
scalar- type gluon propagator D(r) in the coordinate 
space is well described by the Yukawa-type function 
-DYukawa(r) = ime""" '/ V with m = 0.624(8)GeV and 
A = 0.162(2) in the whole range of r = 0.1 ~ l.Ofm. 

As a summary figure, we show in Fig[7]the comparison 
between the obtained Yukawa- type function -DYukawa(r) 
and all the lattice QCD data of the scalar-type gluon 
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FIG. 6: The ratio of the lattice QCD data Aatt(r) at (3 = 
6.0 to the fit-functions -D mass (r), D Yu kawa(»") and Ddipoic(r) 
on the scalar- type gluon propagator, i.e., Diatt /Omasa , 
Aatt/f Yukawa, and -Dlatt/Ai ipolc - 



propagator D(r) at (3=5.7, 5.8, and 6.0 in the range of 
r = 0.1 ~ l.Ofm. For clearer presentation of the Yukawa- 
functional behavior of the lattice gluon propagator D(r), 
we also show the logarithmic plot of rD(r) in Fig[7(b). 
All the lattice QCD data at (3=5.7, 5.8, and 6.0 are found 
to be well reproduced with the Yukawa-type function in 
the range of r = 0.1 ~ l.Ofm. Note also that all the 
lattice data of rD(r) locate around a straight line in the 
logarithmic plot of Fig[7£b). 

Next, we investigate the gluon propagator D(p 2 ) in 
the momentum space in terms of the Yukawa-type func- 
tion. In Fig|51 we show the scalar-type gluon propagator 
D(p 2 ) in the Landau gauge obtained in lattice QCD at 
j3 = 6.0, and L>Yukawa(p 2 ) = 4tt 2 Ato(p 2 +to 2 )" 3 / 2 , which 
is the Fourier transformation of the Yukawa-type func- 
tion -DYukawafV)- The horizontal axis is p = (paPa) 1 ^ 2 - 
Here, we use the same parameters m = 0.624GeV and 
A = 0.162 as those used in the best-fit analysis for the 
coordinate-space gluon propagator. From Figs. and [51 
the lattice QCD data of D(p 2 ) at (3=5.7, 5.8, and 6.0 are 
found to be approximated with I?Yukawa(p 2 ) in the range 
of p < 3GeV. We also perform the best-fit analysis of the 
lattice data of the gluon propagator D(p 2 ) in the momen- 
tum space with I^Yukawa (p 2 ) = 4:Tr 2 Am(p 2 + m 2 )~ 3 / 2 . For 
the lattice data at /?=6.0 in the fit range of p 2 < (3GeV) 2 , 
the best-fit parameters are found to be m ~ 0.577GeV 
and A ~ 0.151, which are close to the values obtained 
from the fit-analysis of the coordinate-space gluon prop- 
agator. 

As a caution, in the UV region of p > 3GeV, the gluon 
propagator D(p 2 ) gradually deviates from the Yukawa 
form -DYukawa^ 2 ) and gradually approaches the pertur- 
bative form 1/p 2 . However, the coordinate-space gluon 
propagator D(r) is found to be almost unchanged for 
r = 0.1 ~ l.Ofm by correcting the UV deviation. 

Also in the Deep-IR region of p < 0.5GeV, as is briefly 
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FIG. 7: The Yukawa-type function Dyukawa (r) = Ame~ mr jr 
(solid line) with m = 0.624GeV and A = 0.162 obtained by 
the fit analysis at /3 = 6.0, and the lattice QCD data of the 
scalar- type gluon propagator D(r) in the Landau gauge at 
13=5.7, 5.8, and 6.0 in the range of r = 0.1 ~ l.Ofm. The 
lower figure is the logarithmic plot of rD(f) and rDy u kawa(r). 



summarized in Appendix B, there appears some devia- 
tion between D(p ) and -Dyukawa [p 2 ) ■ In fact, such a 
Deep-IR deviation is indicated by recent huge- volume lat- 
tice QCD studies [H, |H, [HJ . In the momentum space, 
the true gluon propagator D{p 2 ) turns out to take a sat- 
urated value smaller than the Yukawa-type propagator 
-DYukawa(p 2 ) in the Deep-IR region of p < 0.5GeV. In 
other words, p ~ 0.5GeV is the lower bound on the ap- 
plicability of the Yukawa- type propagator -Dyukawa (p 2 ) 
to the gluon propagator. Based on the huge-volume lat- 
tice data [HJ , we investigate the Deep-IR-corrected gluon 
propagator in the coordinate space in Appendix B. As the 
conclusion, even after the correction in the Deep-IR re- 
gion, the Yukawa-type function is found to work well for 
r = 0.1 - l.Ofm. (See Fig[l4]m Appendix B.) 

As the main conclusion of this paper, we summarize 
the functional form of the gluon propagator in the Lan- 
dau gauge obtained in SU(3) lattice QCD. 

1. The coordinate-space gluon propagator D(r) in 



FIG. 8: The Yukawa-type propagator in the momentum 
space, i.e., Z?Yukawa(p 2 ) = 4n 2 Am(p 2 + m 2 )~ 3 ^ 2 (solid line) 
with m — 0.624GeV and A — 0.162, the same values used in 
FigfT] The horizontal axis is p = (paPa) 1 ^ 2 ■ The symbols de- 
note the lattice QCD data of the scalar-type gluon propagator 
D(p 2 ) in the Landau gauge at (3 = 6.0, where the momentum 
is defined as p M = | sin(^ ii ). 



the Landau gauge is well described by the four- 
dimensional Yukawa-type function as 



D(r) 



_ D aa 

24 w 



(r) 



A — e" 

r 



(42) 



with m ~ 600MeV and A ~ 0.16, for the whole 
region of r = {x a Xa) x l 2 = 0.1 ~ l.Ofm. (This is 
valid even after the possible correction in the UV 
and Deep-IR regions.) 

2. The gluon propagator D(p 2 ) in the momentum 
space is also well described by the correspond- 
ing new-type propagator (four-dimensional Fourier 
transformed Yukawa-type function) as 



Dtp 2 ) = —D aa (p 2 ) = — l ' 
KP ' 24 ""^ 1 (p2- 



2)3/2 



(43) 



with m ~ 600MeV and A ~ 0.16 (same values), in 
the momentum region of 0.5GeV < p < 3GeV. 

Note here that all the component of the gluon propagator 
Df v {x ~y) = (A-(x)A b u (y)} and D%{p) in the Landau 
gauge can be analytically expressed, starting from the 
Yukawa-type function. (See Sec. II.) 

Quantitatively, the Yukawa- type propagator D{r) ex- 
hibits a slower decreasing feature, and D(p 2 ) exhibits 
faster decreasing, in comparison with the ordinary mas- 
sive propagator. It seems suggestive to rewrite D(p 2 ) as 



D(p 2 



Z{p 2 
p 2 + r 



Z(p 2 



4:TT 2 Am 



(P 2 



i 2 y/ 2 ' 



(44) 



where Z(p 2 ) corresponds to the wave-function renormal- 
ization of the gluon field, in a similar manner to the 
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Schwinger-Dyson formalism. Near the on-shell-like con- with r = y/x 2 + 1 2 . Like Eq. ([28l) . the zero-spatial- 



dition of p 2 



= 0, Z(jp 2 ) tends to diverge as 



which leads to anomalous gluon propagation and may 
mimic the gluon confinement, i.e., the absence of on-shell 
gluon states. 

To see the physical meaning of the massive parameter 
to in the Yukawa-type function, we compare the Yukawa- 
type propagator form e~ mr /r with the massive propaga- 
tor form K\(mr)jr in the coordinate space. In spite of 
a significant difference at the short distance, their dif- 
ference is just the prefactor at the large distance, where 
Ki(mr)/r ~ e -™r and the main reduction factor 
is e~ mr . Then, the mass parameter to ~ 600MeV in 
the Yukawa-type gluon propagator directly corresponds 
to the effective gluons mass in the infrared region. (See 
Sec.VI C.) 

Here, we briefly comment on the other functional 
forms for the gluon propagator. Up to now, many func- 
tional forms of the gluon propagator D{p 2 ) have been 
considered and c omp ared with the lattice QCD result 
[H El [H [M [H Hi H, H2 . Some functional forms well 
describe the gluon propagator better than the Yukawa- 
type function, but they need four or more fit parameters 
and take highly non-analytical complicated form. On 
the other hand, the Yukawa- type function £>Yukawa(r) has 
only two parameters {A, to) and takes an analytical form, 
which are its advantages. For example, owing to the an- 
alyticity, the Yukawa- type gluon propagator L>Yukawa(f) 
leads to analytical applications of gluonic nonperturba- 
tive quantities, as will be demonstrated in the next sec- 
tion. 



VI. ANALYTICAL APPLICATIONS OF 
YUKAWA- TYPE GLUON PROPAGATOR 

In this section, as the applications of the Yukawa-type 
gluon propagator, we derive analytical expressions for 
the zero-spatial- momentum propagator Do(t), the effec- 
tive mass M c g(t), and the spectral function p{uj) of the 
gluon field. All the derivations can be analytically per- 
formed, starting from the Yukawa-type gluon propagator 

L>Yukawa(?"). 



A. Zero-spatial-momentum propagator of gluons 

Second, we consider the zero-spatial-momentum prop- 
agator D()(t), associated with the Yukawa- type propaga- 
tor -DYukawa(f)- Do(t) was introduced in Sec.IV-B in the 
context of the effective mass in lattice QCD. Here, we 
mainly deal with the continuum formalism with infinite 
spatial volume. For the simple argument, we first ne- 
glect the temporal periodicity, which is justified for large 
temporal lattice size. 

We start from the Yukawa- type gluon propagator, 



D 



Yukawa 



(r) 



Am 



Am 



V^TT 2 



(45) 



momentum propagator is given by 
1 

24 



D (t) = — I d J x(^(f,i)^(0,0)) = / d'x Z?YukawaM 

(46) 

in the continuum formalism. Using the three-dimensional 
polar coordinate of cc, we calculate this integral as follows: 



D (t) 



An Am I dx x* 
lo 



1 



-my^x 2 -\-t 2 



sfx 2 



t 2 



AjiAin I ilr\ /-'( 



AnAmt" 



dr \J r 2 — le 



= AT:Amt 2 —Ki(mt) = 4irAtKi(mt), (47) 

with f = r/t. Here, we have used Eq. (|A9[) on the mod- 
ified Bessel function. Thus, we derive an analytical ex- 
pression for the zero-spatial-momentum propagator, 



D (t) = ATtAtKximt). 



(48) 



For the actual comparison with the lattice QCD data, 
we take account of the temporal periodicity, which is used 
in lattice calculations. In this case, Do(t) is given as 

D (t) = 4xA[tKi{mt) + (N t - t) Ki(m(N t - t))]. (49) 

In Fig[9l we show the lattice QCD result of Do(t) in 
the Landau gauge, and the theoretical curve of Eq. (|4"9"|) 
with TO=0.624GeV and t4=0.162, the same values used 
in the previous section. The lattice QCD data of Do(t) 
are found to be well described by the theoretical curve, 
associated with the Yukawa-type gluon propagator. 




FIG. 9: The zero-spatial-momentum propagator Do(t) of glu- 
ons in the Landau gauge. The symbols are the lattice QCD 
data at /3 = 6.0, and the solid line denotes the theoretical 
curve of Eq. (|49p . derived from the Yukawa- type propagator 
with m=0.624GeV and A=0.162, the same values in FigE 
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We also consider Do(t) in the lattice formalism, for 
more direct comparison with the lattice QCD data. Here, 
we use the expression of D (t) with -DYukawa(p 2 ), 



D (t) = 





d 4 p 




f°° dpo 


e ipot b- 


Loo 2^ 





"D 



Yukawa 



b 2 ) 



Yukawa 



Cp8) 



(50) 



On the lattice, -Dyukawa(Po) x (Po + ™ 2 ) 3 ^ 2 i s gi yen as 



-DYukawa(Po) 



4ir 2 Ar 



{(2sin(^)) +m 2 } 



3/2 • 



(51) 



with p = 2sin(7rn/iV t ) (n = 0, 1, 2, .... iV f - 1) in the 
lattice unit. Then, we obtain an analytical expression 
for the zero-spatial-momentum propagator, 



, N t -1 

Do® = Jf t V 



At: 2 Am 



n=0 



{(2sin(^)) 2 + m2} 



3/2 ' 



(52) 



No significant numerical difference turns out to be found 
between the continuum and the lattice formulae, under 
the condition of our lattice QCD calculation. This means 
that the integral quantity Do(t) is not so sensitive to the 
details of the UV-behavior of the propagator -DyukawaM- 
In the UV region, the Yukawa-type propagator 
-DYukawa(p 2 ) °c (p 2 + m 2 ) -3 ' 2 deviates from the cor- 
rect behavior of the perturbative propagator 1/p 2 . For 
the quantitative estimate of the influence from the devi- 
ation in the UV region, we calculate the zero-spatial- 
momentum propagator Do(t) using the UV-corrected 
Yukawa propagator D^^l(p 2 ), which is D Yukawa (p 2 ) 
for p < 4GeV and 1/p 2 in the UV region of p > 4GeV. 
The difference of Do(t) is found to be at most 1% between 
the cases with ^Yukawab 2 ) and D^Z(p 2 )- Thus, the 
integral quantity D (t) is insensitive to the UV behavior 
of the gluon propagator, so that it is also expected to be 
insensitive to the discritization error [25|, [2g, [54| in the 
UV region in the lattice calculation. 



B. Effective mass of gluons 

Second, we investigate the effective mass M c g(t), as the 
consequence of the Yukawa- type propagator -DyukawaM- 
The effective-mass plot is a general useful technique for 
the mass estimation in lattice QCD, and was actually 
examined for gluons in Sec.IV-B. For simplicity, we here 
treat the three-dimensional space as a continuous infinite- 
volume space, while the temporal variable t is discrete 
and periodic. For the simple notation, we here use the 
lattice unit for t. 

When the temporal periodicity can be neglected, the 
zero-spatial- momentum propagator Do(t) is expressed by 
Eq.(|48p. and we obtain an analytical expression of the 



effective mass, 



M cff (t)=ln T ^ TT =l::, 



tK x {mt) 



D (t + 1) {t+\)K x {m{t+\)) 



(53) 



From the asymptotic form of K\{z) cx z 1 / 2 e z , the ef- 
fective mass of gluons is approximated as 



Meff(i) si m - i hi (1 + i) 



for large t. This functional form indicates that M e g(t) is 
an increasing function and approaches m from below, as 
t increases. 

Note that the mass parameter m in the Yukawa-type 
gluon propagator directly corresponds to the effective 
mass M e g(t) of gluons for large t. In fact, m ~ 600MeV 
has a definite physical meaning of the effective gluon mass 
in the infrared region. 

Note also that the simple analytical expression of 
Eq. (|53p or (|54|) reproduces the anomalous increasing be- 
havior of the effective mass M e g(t) of gluons, as was ob- 
served in FigJU As for the increasing behavior, there is a 
general argument: this can occur when the spectral func- 
tion is not positive definite [TH, HO, H3, [H[ , although the 
concrete form of the gluon spectral function is not yet 
known. Instead, this framework with the Yukawa-type 
gluon propagator gives an analytical and quantitative 
method, and is found to well reproduce the lattice result. 
(The actual comparison is demonstrated with Eq. (f55]) .) 

Next, we take account of the temporal periodicity, 
which is used in lattice QCD calculations. In this 
case, the effective mass M e g(t) is defined by "cosh- type" 
as Eq. (l32j) . and the zero-spatial-momentum propagator 
D (t) is given by Eq. (|4*§|) . Then, the effective mass 
M c tf(t) of gluons is expressed as 

cosh[M eff (t)(iV t /2-(t + l))] 



(55) 



cosh[M eS (t)(N t /2-t)} 
(t + l)gi(m(j + 1)) + (N t -(t + lJJK^mjNt - (t + 1))) 
tK^mt) + (N t - t)ifi(n»(Wi - t)) 

In FigfTOl we show the theoretical curve obtained by 
Eq. (|55]) together with the lattice result of M e g(t). Here, 
we take m=0.624GeV, the same value used in the previ- 
ous section. The lattice QCD data of M e gp(t) are found 
to be well described by the theoretical curve, associated 
with the Yukawa-type gluon propagator. 



C. Spectral function of gluons in Landau gauge 

From the analytical expression of the zero-spatial- 
momentum propagator Do(t) = AmAtK\(m£) in Eq. (|4"5)l . 
we can derive the spectral function p(io) of the gluon 
field, associated with the Yukawa-type gluon propaga- 
tor. For simplicity, we take continuum formalism with 
infinite space-time. 
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FIG. 10: The effective mass M e s(t) of gluons in the Landau 
gauge. The symbols denote the lattice QCD data at f3 = 6.0, 
and the solid line denotes the theoretical curve of Eq. (|55[) 
derived from the Yukawa-type propagator with m=0.624GeV, 
the same value used in Fig[7] 



The relation between the spectral function p(co) and 
the zero-spatial- momentum propagator Dq (t) is given by 
the Laplace transformation, 



Do{t) = I duj p(ui) e 
Jo 



(56) 



When the spectral function is given by a ^-function such 
as p{uj) ~ S(u> — wo), which corresponds to a single mass 
spectrum, one finds a familiar relation of Do(t) ~ e - " '. 
For the physical state, the spectral function p{to) gives 
a probability factor, and is non-negative definite in the 
whole region of to. This property is related to the unitar- 
ity of the S-matrix. 

In general, the Laplace transformation is defined by 



g(t) = / du; e~ M f(u), 
Jo 



(57) 



and the inverse Laplace transformation is expressed as 

rc-\-ioc 



(58) 



Then, from Eq. (fS"6")) , the spectral function p(oj) is ex- 
pressed as 



p(uj) = — dt e ut D (t) 

ZTTt J c—ioo 
-i pc-\-ioo 

— I dt e ut ATtAtKx{mt) 
27ri J c—ioo 



2m 
1 

2TTi 



c—too 

C ' +lC ° V 4ttA . . . . 

<:~* — -t'Ki(t'), (59) 



dt' ^ * 



C — ZOC 



111 ' 



with a/ = uj/m, t' = mt and c' — mc. 



Performing the partial integration in Eq. (|A9D on the 
modified Bessel function, we obtain a formula of the 
Laplace transformation, 



Ki(t) 



doj e 



dto 



to 



(iU 2 - 1)V2 



W -l), (60) 



(^-1)1/2 

which leads to the inverse Laplace transformation, 

i pc+ioc 

^ii^ dte "" Kllt) = W^W^- l) ' (61) 

By differentiating this formula by we find 

^ pc-\-ioo 



2m 



dt e wt 



tK^t) 



c—too 

1 



LO 



= -(^W^- 1)+ w=w* * {u - 11 (62) 

where the second term includes an infinite factor besides 
the (^-function. Then, we apply this formula to Eq. (|59[) . 
and obtain the spectral function p(uj) as 



_L / dt'e^ ^t'K^t'l 



- - 7 ) + ^^S^-m). (63) 

For more rigorous derivation, we avoid the singularity 
at u> — m by regularizing Eq. (|6ip as 

-J c ioo d te^K l{ t) = w - w B { .-l-e) (64) 



with a positive infinitesimal £, and find the formula of 

pc-\-ioo 

2m 



-i pc-\-ioo 

— dt e ut tKi(t) 

Jin Jc—iso 



1 



-6(lo 



1 



( w a _ 1)3/2 

which leads to the regularized spectral function 
An Am 



(65) 



p e (w) = - 



4ttA 

— ——^u-rn-e ) + — S{uj-m-e). 

[lu 1 — m z yi z (2me) L >' ! 

(66) 

In the calculation of the Laplace transformation from 
p(u>) to D{j{t), we can avoid the divergence of the inte- 
gral at to = m by using the regularized spectral function 
p £ (w). Then, we can perform the integration, and prop- 
erly obtain Do(t) = 4TrAtKi(mt), by taking the limit of 
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e — > after the integration. We have numerically con- 
firmed that the spectral function p s (u>) in Eq. (|66p repro- 
duces Do(t) = 4:irAtKi(mt) by the Laplace transforma- 
tion in the limit of e — > 0. 

In this way, we derive the spectral function p(u>) of the 
gluon field, associated with the Yukawa-type propagator: 



p(w) 



4ttAt 



iff 



2)3/: 



r 0(w - m) + n72 (5( - w ~ m )' 

(w — mj 1 '^ 

(67) 

which is regularized as p £ (uj) in more rigorous derivation. 
Here, m ~ 600MeV is the mass parameter in the Yukawa- 
type function for the Landau-gauge gluon propagator. 
The first term expresses a negative continuum spectrum, 
and the second term a ^-functional peak with the residue 
including an infinite factor, which is positive as e -1 / 2 at 
to = m + £. 

We show in FigJTT] the spectral function p(uS) of the 
gluon field. Although the appearance of the negative- 
value region in the gluon spectral function is expected, 
p(u>) exhibits two anomalous behaviors: it has a positive 
S- functional peak with the residue of +00 at to = m (+e), 
and it takes negative values for all the region of lu > m. 
This negative contribution of the spectral function gives 
unusual behavior of the effective mass M e g(t) of gluons, 
i.e., its increasing behavior on t. 



p(co) 



CO [m] 

FIG. 11: The spectral function p(ui) of the gluon field, as- 
sociated with the Yukawa-type propagator. The unit is nor- 
malized by the mass parameter m ~ 600MeV. As Eq. (|67[) 
indicates, p(u)) shows anomalous behaviors: it has a positive 
8- functional peak with the residue of +00 at u) = m (+e), and 
takes negative values for all the region of u) > m. 



As was discussed in Sec. IV, if the spectral function is 
non-negative definite, like the ordinary hadronic correla- 
tor, the effective mass must be a monotonously decreas- 
ing function, due to the mathematical nature of the sum- 
mation of e~ ut with non-negative coefficients, as can be 



proven with Eq. (|34|) . Physically, this is due to a larger 
positive contribution of the excited states to the effec- 
tive mass at a shorter distance. However, if the spectral 



function includes the ne gati ve- value region, such a defi- 
nite tendency is lost [H S3, 111 EHH • This is the case 
of the gluon field in the Landau gauge. 

More precisely, from Eq. (|54"]) , we find that the increas- 
ing property of M e g(t) can be realized, only when there 
is some suitable coexistence of positive- and negative- 
value regions in the spectral function piui). As a remark- 
able fact, the obtained gluon spectral function p(u>) is 
negative-definite for all the region of lu > m, except for 
the positive ^-functional peak at u) = m. The negative 
property of the spectral function in coexistence with the 
positive peak leads to the anomalous increasing behavior 
of the effective mass M e s(i). 

In the UV region, the true gluon propagator deviates 
from the Yukawa-type function. Then, we estimate the 
contribution from the UV part of this spectral function 
p(lo) to the zero-spatial- momentum propagator D${t). 
Introducing the UV cutoff A, we define the integration 



p(uj) 



(68) 



which reduces to the Laplace transformation in the A — > 
00 limit. For t > O.lfm, the A-dependence of the integral 
quantity Do(t) is found to be negligible for A > 2 GeV, 
i.e., Dq (t) ~ D (t). In other words, due to this UV- 
insensitivity, the obtained information on the UV part 
of the spectral function suffers from a large uncertainty, 
as is also seen in the maximum entropy method (MEM) 
analysis. On the other hand, the IR/IM part of the spec- 
tral function p(u>) is relatively stable and reliable. 

We note that the gluon spectral function p{u) is diver- 
gent at lu = m + £, and the divergence structure is com- 
plicated and consists of two ingredients: a ^-functional 
peak with a positive infinite residue and a negative wider 
power-damping peak. On the finite-volume lattice, these 
singularities are to be smeared, and p(uj) is expected to 
take a finite value everywhere on to. On the lattice, we 
conjecture that the spectral function p{to) includes a nar- 
row positive peak stemming from the 5-function in the 
vicinity of lu = m (+e) and a wider negative peak near 
lu ~ m in the region of lu > m. 

In this way, the Yukawa-type gluon propagator sug- 
gests an extremely anomalous spectral function of the 
gluon field in the Landau gauge. Note that this frame- 
work gives an analytical and concrete expression for the 
gluon spectral function p(co) at the quantitative level. 
Actually, the resulting effective mass M e e(t) well de- 
scribes the lattice result, as shown in FigJU The obtained 
gluon spectral function p(u>) is negative almost every- 
where, and includes a complicated divergence structure 
near the "anomalous threshold", lu — m (+e). 

These anomalous features of the gluon spectral func- 
tion may have some relation to the various nonperturba- 
tive QCD phenomena, such as the gluon confinement and 
the gluonic instability of the QCD vacuum, e.g., gluon 
condensation, the Savvidy vacuum [55[, and the Copen- 
hagen vacuum [56]. In any case, the Yukawa-type gluon 
propagator -DyukawaM, which models the Landau-gauge 
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gluon propagator, would be useful for the analytical and 
quantitative investigation of nonperturbative QCD. 



VII. SUMMARY AND DISCUSSIONS 

We have studied the gluon propagator D^ v (x) in the 
Landau gauge in the infrared/intermediate region of r = 
(x^x^) 1 / 2 = 0.1 ~ l.Ofm, which is relevant to the quark- 
hadron physics, in SU(3) lattice QCD at j3 = 5.7, 5.8, and 
6.0 at the quenched level. From the gluon propagator 
analysis and the effective-mass plot, we have estimated 
the effective gluon mass of 400 ~ 600MeV in the infrared 
region of r = 0.5 ~ l.Ofm. The effective gluon mass 
exhibits a significant r-dependence in this region: it takes 
a smaller value in the smaller r-region. 

We have also studied the functional form of the gluon 
propagator D"^(x) in lattice QCD. As a remarkable fact, 
the lattice QCD result of the Landau-gauge gluon prop- 
agator D^(r) is fairly well described by the Yukawa- 
type form Z?Yukawa(f) oc e~ mr /r with the mass parame- 
ter to ~ 600MeV in the whole region of r = 0.1 ~ l.Ofm 
in four-dimensional Euclidean space-time. This Yukawa- 
type propagator corresponds to the new-type propagator 
£>Yukawa(p 2 ) oc (p 2 + m 2 )~ 3 / 2 in the momentum space, 
through the Fourier transformation, and this also well 
describes the lattice QCD result of the gluon propagator 
Dftftip 2 ) in the momentum space. 

As the application of the Yukawa-type gluon propa- 
gator, we have derived the analytical expression of the 
zero-spatial- momentum propagator Dq (t) , and the effec- 
tive mass M e ff(t). The obtained analytical functions for 
D (t) and M e g(t) well reproduce the lattice QCD re- 
sults, in particular the anomalous increasing behavior of 
M e g(t). We have found that the mass parameter to of 
the Yukawa-type gluon propagator directly corresponds 
to the effective gluon mass in the infrared region of ~lfm. 

We have also derived the analytical expression of the 
spectral function p(ui) of the gluon field, associated 
with the Yukawa-type gluon propagator, using the in- 
verse Laplace transformation of the temporal propaga- 
tor Do(t). As a remarkable fact, the obtained spectral 
function p(w) is negative-definite almost everywhere for 
lo > to, except for a positive ^-functional peak with the 
residue of +oo at lo = to. The coexistence of negative- 
and positive-value regions of p(to) lead to the anomalous 
increasing behavior of the effective mass M s g(t) of glu- 
ons. Thus, the theoretical analysis with the Yukawa-type 
gluon propagator gives a new analytical and quantitative 
method for the nonperturbative gluonic phenomena. 

The Yukawa function -DyukawaM oc e~ rnr /r in the co- 
ordinate space is highly analytic, and this analyticity 
plays an important role in deriving the analytical ex- 
pressions for Do(t), M e s(t), and p(u>). On the other 
hand, the momentum-space propagator -Dy u kawa(p 2 ) oc 
(p 2 +to 2 )~ 3 / 2 includes a singular cut stemming from the 
square root. If this singularity is taken seriously, there 
would arise a problem, because analyticity of the Green 



function is important in quantum field theories, e.g., 
in the Wick rotation converting between the Euclidean 
space and the Minkowski space. Of course, this Yukawa- 
type function is an approximate form for the Landau- 
gauge gluon propagator in the region of r = 0.1 ~ lfm. 
There would be more regular and more thorough expres- 
sion for the gluon propagator. 

In this paper, we have mainly considered the gluon 
propagator in the coordinate space instead of the momen- 
tum space, since the coordinate-space variable is more di- 
rectly obtained in lattice QCD. For the confined particles, 
however, the reason to use the momentum representation 
would be less clear, compared with ordinary particles. In 
the ordinary particles, the momentum representation of 
the Green function is clearly useful to express the pole 
structure, to distinguish the on-shell and off-shell states, 
and so on. However, for the confined field, there is no 
on-shell state, i.e., no physical asymptotic state, so that 
there is no definite reason to use the momentum rep- 
resentation, besides the total momentum conservation. 
Indeed, it is difficult to image the non-zero momentum- 
space propagator with no pole, and the coordinate-space 
representation may be more convenient for some descrip- 
tion of the confined field, similar to potential problems 
in quantum mechanics, where the coordinate-space wave- 
function is convenient. 

The Yukawa- type gluon propagator -DyukawaM in- 
cludes e~ mr as the main reduction factor in the infrared 
region, and m ~ 0.6GeV can be regarded as the infrared 
effective gluon mass. In terms of the infrared reduction, 
a simple "constituent gluon picture" may be approxi- 
mately obtained as Mqb — 2m for the glueball mass 
Mqb- In general, the glueball mass Mqb can be esti- 
mated from the infrared reduction of the glueball corre- 
lator Ggb{x — y) = (3>gb(z)$gb(z/)) with the glueball 
operator, e.g., $gb = G" G" for the scalar glueball, 
with the field strength tensor G M „ [1, [H|. By the Wick 
contraction, Gqb{%) can be expressed by some combi- 
nation of the gluon Green functions. In the framework 
with the Yukawa-type gluon propagator, the leading re- 
duction term of Gqb(^) can be expressed with some 
derivative of {-Dy u kawa(?")} 2 for large r. Then, Gqb(x) 
includes e~ 2mr as one of the main infrared reduction fac- 
tors, and the lowest-glueball mass is roughly estimated 
as Mqb — 2to = 1.2 ~ 1.3GeV, by neglecting all the 
prefactor and higher-order terms, which express the in- 
teractions between gluons. In spite of the crude estimate, 
this value gives the same order of the lowest-glueball mass 
of about 1.5GeV obtained in lattice QCD [l,[44j]. 

In this subject, there remain difficult problems related 
to the confinement mechanism, the infinite- volume limit, 
the Gribov copies, and the gauge dependence. For ex- 
ample, it is a highly difficult mathematical problem to 
find out the precise description of the physical hadronic 
states like glueballs in terms of the confined particles, 
quarks and gluons. As other example, Zwanziger's the- 
orem jl5[ is derived from the argument of the Gribov 
horizon: all connected gluon correlation functions includ- 
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ing the gluon propagator must vanish at zero momentum 
in the infinite-volume limit. However, our lattice QCD 
results and the Yukawa- type propagator D(p 2 ) indicate 
an infrared non-vanishing property of the gluon propa- 
gator as D(p 2 = 0) = 4ir 2 A/m 2 . This may be due to 
the absence of an infinite-volume effect in the Deep-IR 
region, as is conjectured by analytical studies [TEI, Il7j |. 
For this problem, it is desired to clarify the Deep-IR be- 
havior of the gluon propagator, and recent huge- volume 
lattice studies [H, [H, [13] and a recent analytical study 
based on the Schwinge-Dyson equation [l!| also indicate 
the infrared non- vanishing property of the gluon propaga- 
tor. As for the Gribov-copy problem, it is reported that 
the Gribov-copy effect is quantitatively rather small in 
the actual lattice QCD calculation for the Landau-gauge 
gluon propagator [24], [2{|, while the ghost propagator 
slightly suffers from it. However, this is a fundamental 
problem in QCD, and it would be serious in the argument 
of the large- volume limit, so that it is also desired to re- 
move the Gribov copies and to extract the fundamental 
moduli region. 

Finally, we discuss the Yukawa-type gluon propagation 
and a possible dimensional reduction due to the stochas- 
tic behavior of the gluon field in the infrared region. As 
shown in this paper, the Landau-gauge gluon propaga- 
tor is well described by the Yukawa function in four- 
dimensional Euclidean space-time. However, the Yukawa 
function e~ mr /r is a natural form in i/iree-dimensional 
Euclidean space-time, since it is obtained by the three- 
dimensional Fourier transformation of the ordinary mas- 
sive propagator (p 2 + to 2 ) -1 . In fact, the Yukawa- type 
propagator has a "three-dimensional" property. In this 
sense, as an interesting possibility, we propose to inter- 
pret this Yukawa-type behavior of the gluon propagation 
as an "effective reduction of the space-time dimension" . 

Such a "dimensional reduction" sometimes occurs in 
stochastic systems, as Parisi and Sourlas pointed out for 
the spin system in a random magnetic field On the 
infrared dominant diagrams, the D-dimensional system 
coupled to the Gaussian-random external field is equiv- 
alent to the (D — 2)-dimensional system without the ex- 
ternal field. In fact, the space-time dimension of the the- 
ory is apparently reduced by two. For the system cou- 
pled to the Gaussian-random external source, the dimen- 
sional reduction is universal and is associated with a hid- 
den supersymmetry: in the superspace formalism with 
(x^,6, 9), the integration over two Grassmann- variables 
(6,6) reduces the space-time coordinates x M by two 

We note that the gluon propagation in the QCD vac- 
uum resembles the situation of the system coupled to the 
stochastic external field. In fact, as is indicated by a large 
positive value of the gluon condensate (G a ^ u G a ^ u ) > in 
the Minkowski space, the QCD vacuum is filled with a 
strong color- magnetic field [f| H, HU, [56| , which can con- 
tribute spontaneous chiral-symmetry breaking (58j . and 
the color-magnetic field is considered to be highly random 
at an infrared scale [56|, H^, [6(3|. Since gluons interact 
each other, the propagating gluon is violently scattered 



by the other gluon fields which are randomly condensed 
in the QCD vacuum at the infrared scale, as schemati- 
cally shown in Figfl2l 




ng gluon 



FIG. 12: A schematic figure for a propagating gluon in the 
QCD vacuum. The QCD vacuum is filled with color-magnetic 
fields which are stochastic at an infrared scale, and the gluon 
propagates in the random color-magnetic fields. 

Actually in the infrared region, the gluon field shows 
a strong randomness due to the strong interaction, and 
this infrared strong randomness is considered to be re- 
sponsible for color confinement, as is indicated in strong- 
coupling lattice QCD @, Q . In the chiral random matrix 
theory for QCD [59|], the infrared randomness of gluons is 
taken account as an essence of QCD in a simplified man- 
ner, and the gluon field is replaced by a Gaussian-random 
external field coupled to the quark field. 

Notice that there are two kinds of randomness in the 
gluon field: one is a completely-random gauge degree of 
freedom, which is fake, and the other is a net physical 
gluonic degree of freedom, which is not completely but 
highly random at the infrared scale. For the argument 
of physical randomness, these two concepts should be 
clearly divided, since the gauge degree of freedom is just 
a fake. Even after the removal of the fake gauge degree of 
freedom by g aug e fixing, the gluon field exhibits a strong 
randomness [601 ] accompanying a quite large fluctuation 
at the infrared scale. 

As a generalization of the Parisi-Sourlas mechanism, 
we conjecture that the infrared structure of a theory in 
the presence of the quasi-random external field in higher- 
dimensional space-time has a similarity to the theory 
without the external field in lower-dimensional space- 
time. From this point of view, the Yukawa-type behavior 
of gluon propagation may indicate an "effective reduction 
of the space-time dimension" by one, due to the stochas- 
tic interaction between the propagating gluon and the 
other gluon fields in the QCD vacuum, of which net phys- 
ical fluctuation is highly random at the infrared scale. 

In this paper, the Yukawa-type gluon propagator is 
obtained phenomenologically from lattice QCD results, 
but we expect some deeper theoretical reasons for the 
Yukawa-type propagation, which may be an effective re- 
duction of the space-time dimension, due to the stochas- 
tic behavior of the infrared gluon field. In any case, the 
Yukawa-type gluon propagator would provide a new ana- 
lytical framework for the study of nonperturbative QCD. 
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APPENDIX A: FOURIER TRANSFORMATIONS 



and its inverse Fourier transformation gives Eq. (|35[) . 



2. The massive propagator in four- and 
two-dimensional Euclidean space-time 

We here calculate the Fourier transformation of the 
massive propagator (p 2 + m ) . 



In this Appendix, we derive several Fourier transfor- 
mations used in SecfVl 



1. The Yukawa function in the four-dimensional 
Euclidean space-time 

We first derive the four-dimensional Fourier transfer- 
or by calculating 

p—mr 

/Yukawa = / #X e"*- , (Al) 



with 



(x^x^) I . 



We use the polar coordinate 



(r, 9q , 6i, 82) polar m four-dimensional Euclidean space- 
time, and choose the axis to satisfy p ■ x = pr cos 82, 
without loss of generality. Then, the Fourier integral 
(|Aip is expressed as 



'Yukaua = 4-7T / dd 2 SIB. 9 2 

>0 JO 

= Air I d8 2 sin 2 8 2 



dr r 2 e -( m - i P cos9 2) r 
1 



8tt / d8 2 sin 2 8 2 - 



(m — ip cos 9 2 ) 3 J 
1 

m — ipcos^) 3 



dt i 2 e~* 



(A2) 



Here, we have replaced (m — ip cos 9 2 )r by t, and changed 
the integration range of t using the analytic continuation. 
We rewrite /yukawa with partial integration as 



* Yukawa 



/ dusmv— 

p J Q d9 

4iri 



1 



(m — ip cos 9) 2 
1 



d9cos9-— . 
p J (m — ipcos8) z 

1 



P dp J m — ip cos ( 

47T d 7T 



4tt 2 



p dp (p 2 + to 2 ) 1 / 2 (p 2 + m 2 ) 3 / 2 
Here, we have used the integral formula, 

1 7T 



(A3) 



d8- 



a + ibcos8 y/a 2 



b 2 



(a,beR). (A4) 



Thus, we obtain the Fourier transformation of the 
Yukawa function as 



^Yukawa = / d i Xe W ' X - 



4tt 2 



(p 2 + m 2 ) 3 / 2 



(A5) 



a. Four-dimensional Euclidean space-time 

First, we consider the Fourier integral of (p 2 + m 2 ) -1 
in the four-dimensional Euclidean space-time, 



'4dim — 



d4 P jp-x L 



(2tt) 4 p 2 + m 2 



(A6) 



By rotating the coordinate, we set x — (r, 0, 0, 0) without 
loss of generality. Then, the integral is expressed as 



(2tt) 3 J_ x 2ir~ pl+jP + m 2 
(2tt) 3 2y / p 2 +m 2 



(A7) 



with p = (pa , p) . With the three-dimensional polar coor- 
dinate of p, /4dim is written as 



^4di: 



dpp 



2 1 \ ,"--! m-r 



4ir 2 J " y/ p 2 + m 
1 f°° 

-I dE^ 2 

n J m 



-y/p 2 +r. 



4?T 



m 2 e- Er 



— 5-m 2 / dt\J e 2 — It 
4?H J 1 



(A8) 



with E = \J p 2 + m 2 and e = E/m. Using the integral 
representation of the modified Bessel function, 

/>00 

K 1 {z)=z dte- zt (t 2 - 1) 1/2 (Re z > 0), (A9) 



we obtain the Fourier transformation formula, 



'4dim 



dA V i P . x 1 



1 TO 



(27r) 4 p 2 + to 2 47T 2 r 



= T~2 — Ki(mr). (A10) 



b. Two-dimensional Euclidean space-time 

Next, we consider the Fourier integral of (p 2 + m 2 ) -1 
in the two-dimensional Euclidean space-time, 



1 



d 2 p 

(27r) 2 p 2 + TO 2 



(All) 
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By rotating the coordinate, we set x = (r, 0) without loss 
of generality, and integrate P2 as 



dpi 

dpi < 



dp 2 



1 



x> 271 J-oo 2n 
ipi r 



Pl+P2+ m2 



dpi cos(pir) 



,(A12) 



Using Mehler's integral representation of the modified 
Bessel function, 



. . r°° , cos(zt) 

we obtain the Fourier transformation formula, 



(A13) 



d 2 p 



p 2 + m 2ir 



K (mr). (A14) 



3. The dipole-type propagator 

We deal with the Fourier integral of the dipole-type 
propagator (p 2 +m 2 )~ 2 in the four-dimensional Euclidean 
space-time, 



id! 



dipolc 



1 



d 4 v 

(2tt) 4 \p 2 + m 2 ) 2 



1 d 

2m dm 



d A p 



1 



(2tt) 4 p 2 + m 2 
Using the Fourier transformation (|A10|) . we get 

1 m 



(A15) 



-^dipolc 



1 d 
2m dm 
1 

"8^ 



4.71-2 r 



K\ (mr) 



K\(mr) H — K\(mr) 

mr r dm 



(A16) 



From the relation of the modified Bessel function, 

zK' v {z) + vK v (z) = -zK v -i{z), (A17) 
we obtain the Fourier transformation formula, 



i, 



dipolc 



d 4 P 

(27T) 4 



1 



1 



(p 2 + m 2 ) 2 Sir 



2 K (mr). 



(A18) 



APPENDIX B: DEEP-IR CORRECTED GLUON 
PROPAGATOR 



coordinate-space gluon propagator D(r) in the IR/IM 
region of r = 0.1 ~ l.Ofm. 

Figure [13] shows the scalar-type gluon propagator 
D(p 2 ) in the recent lattice-QCD calculation with a huge 
volume 96 4 at (3 — 5.7, taken from Ref.[32|. Here, 
a renormalization constant is multiplied for the huge- 
volume lattice data so as to adjust them to the renor- 
malization condition (|20[) at /i=4GeV. In the momentum 
space, the true gluon propagator D(p 2 ) turns out to take 
a saturated value smaller than the Yukawa-type propa- 
gator I?Yukawa(p 2 ) in the Deep-IR region of p < 0.5GeV. 
In other words, p ~ 0.5GcV is the lower bound on the 
applicability of the Yukawa- type propagator £>Yukawa(p 2 ) 
to the gluon propagator. 
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FIG. 13: The infrared behavior of the gluon propagator 
D(p 2 ). The triangle symbols denote recent huge- volume lat- 
tice data taken from Ref. '32]. The solid line denotes the 
Yukawa- type propagator .DYukawa(p 2 ), and the dashed line 
the Deep-IR-corrected propagator -DYukawa(p 2 ) with pm. = 
0.45GeV. 



Taking account of the deviation in the Deep-IR region, 
we define the Deep-IR-corrected momentum-space prop- 
agator as 



r) !l;.,,n _\ I D Yukawa (p ! /' PlB . m , 

[-DvukawabiRj (COnSt.) V < PlR 



with the IR-saturation momentum of pir = 0.45GeV. 
This value of pir is determined so as to consist with the 
huge-volume lattice result in the Deep-IR region. Us- 
ing this Deep-IR-corrected propagator ^Yukawa (Pir)j we 
calculate the Deep-IR-corrected coordinate-space propa- 
gator by the Fourier transformation as 



In the Deep-IR region, there has been reported to 
appear some deviation on the momentum-space gluon 
propagator D(p 2 ) between small-size lattice data and 
huge- volume lattice data [H, [H, HH • In this Appendix, 
we demonstrate that the deviation in the Deep-IR re- 
gion does not affect the Yukawa-type behavior of the 



nlRcorr / 
^Yukawa V 



(2ny 



-ip-X j~y 



IRcorr 
Yukawa 



(P 2 



(B2) 



In Figrni we show the Yukawa-type propaga- 
tor -DYukawaM an d this Deep-IR-corrected propagator 
D^ a (r), together with the lattice QCD data. The 
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FIG. 14: The Yukawa- type propagator DYukawa(r) (solid line), 
and Deep-IR-corrected propagator I^Yukawa ( r ) (dashed-line) , 
together with the lattice data. The difference between them 
is fairly small in the IR/IM region of r = 0.1 ~ l.Ofm. 



Thus, in the coordinate space, the Yukawa-type func- 
tion I?Yukawa(^) works well for the IR/IM region of r = 
0.1 ~ 1.0 fm, even after the correction in the Deep-IR 
region. 
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